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Abstract 

The poor angular resolution of imaging 7-ray telescopes is offset by the large re- 
flector areas of next generation telescopes such as MAGIC (17 m diameter), which 
makes the study of optical emission associated with some 7-ray sources feasible. 
Furthermore, the extremely fast time response of photomultipliers (PMs) makes 
them ideal detectors for fast (subsecond) optical transients and periodic sources 
like pulsars. The optical pulse of the Crab pulsar was detected with the HEGRA 
CT1 central pixel using a modified PM, similar to the future MAGIC camera PMs. 

The HEGRA CT1 telescope is a 2.1 m radius Cerenkov telescope, which has 
seen the Crab optical pulsations. The corresponding observation time required for a 
detection is 200 seconds, but MAGIC will only require ~ 30 sec if the same strategy 
with standard PMs is used. 

The constancy of the pulsed fraction (with a relatively small error) during the 
2 hours CT1 measurements shows that the pointing did not change during the 
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observations. The purpose of these periodic observations was to determine the light 
of the night sky (LONS) for the galactic anticenter Crab region, and to compare 
it with the LONS for a nearby bright star (£ Tauri). These obtained LONS values 
vary between 4.6 and 4.3 xlO 12 ph ■ m~ 2 • s" 1 • sr" 1 (with a coarse estimated error 
of 10%), being the first value from the timing signature of the Crab pulsar, and 
the second one was derived from the technique using the ( Tauri standard star. 
Our results are between 2.5 and 3 times larger than the previously measured LONS 
(outside the galactic plane) of (1.7 ± 0.4) x W 12 ph ■ mT 2 ■ s _1 • sr _1 , as expected 
since the Crab pulsar is in the galactic plane, which implies a slightly higher energy 
threshold for Crab observations, if the higher value of CT1 measured LONS rate 
for galactic sources is used. 
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1 Introduction 

Imaging Atmospheric Cerenkov Telescopes (IACTs) can be used to detect the 
optical emission of an astronomical object through the increased DC current 
of the camera pixels. Each of these pixels has a field-of-view (FOV) in the 
range of 0.1-0.3°. Hence, when an IACT follows a star, its optical emission is 
fully contained in the telescope central pixel. 

The Crab Nebula is known to be the most important calibration source for 

Email addresses: emma@gae.ucm.es (E. Oha-Wilhelmi) , cortina@ifae.es (J. 
Cortina), fskocdj@puk.ac.za (O.C. de Jager), fonseca@gae.ucm.es (V. 

Fonseca) . 



2 



high energy astrophysics, and especially for the ground-based atmospheric 
IACTs. The night sky background (LONS) influences the trigger threshold, the 
energy and the spectral determination of 7-rays, therefore the LONS around 
the Crab Nebula must be accurately known in the wavelength interval of the 
photomultipliers (PMs) used for the IACT cameras. Mirzoyan and Lorenz J?| 
measured the LONS at the site of the HEGRA experiment for a narrow angle, 
using a detector attached to the telescope CT2, which integrates only the 
starlight from the background of faint stars, obtaining a value of 7iOA r s = (L7± 
0.4) x 10 12 ph ■ mT 2 ■ s^ 1 ■ sr^ 1 (outside the galactic plane). In this paper we 
measure this Ilons for the Crab region and we compare it with the value of 
J?|), using a modified PM in the central pixel of the CT1 camera. 

The optical pulsed emission of the Crab pulsar has already been detected using 
IACTs (P) and other Cerenkov detectors (4). The main purpose of such studies 
are to determine the timing parameters of the pulsar, with which the 7-ray 
arrival time can be folded, thus reducing the search for pulsed emission to a 
single statistical trial. It also allows us to test the accuracy of the used timing 
system and the solar system barycentric correction software, when comparing 
the timing parameters against contemporary radio data. This paper reports 
the optical timing signature of the Crab pulsar obtained, using the central pixel 
of CT1 and the DC observations of the bright nearby star ( Tauri, used to 
derive the LONS around the Crab Nebula. These two values cross-check the 
LONS measurements obtained using two different techniques: temporal (Crab 
pulsation) versus spatial (nearby star). 

Our results with CT1 can be extended to a much larger telescope such as 
MAGIC (1) (17 meter diameter). IACTs are calibrated using as reference the 
Crab signal, and we can therefore use this regular Crab monitoring, not only 
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for Cerenkov measurements, but also to measure the optical Crab pulsation 
with a modified PM in the telescope camera. Furthermore, since the astro- 
nomical community does not monitor the Crab pulsar continuously (on a 24 
hour timescale) in radio or optical, observations with Cerenkov telescopes will 
also help to serendipitously discover pulsar glitches or period irregularities. 

2 The Crab Optical Spectrum and Expected Single Photoelectron 
Response. 

The optical spectra of the Crab pulsar main pulse and interpulse, and the 
spectrum of the underlying nebula (J2I), follow the power laws given in table 1. 
The nebula with a size of 4' x 6' fits well into a single pixel. 



Flux 


K a 


Main pulse 
Interpulse 
Underlying nebula 


5.9 X 10" 15 +0.2 ± 0.1 
1.9 x 10" 15 +0.2 ± 0.1 
8 x 10~ 12 -0.4 ± 0.1 



Table 1 

Flux in the 0.1 phase window centered on the main pulse and interpulse, where, 
F X =K x (X/X )- a - 2 , X = 6000 A, 5000 A < A < 7500 A. The useful wavelength 
range is 3000-8000 A. Fluxes are in erg ■ s" 1 • A" 1 ■ cm~ 2 . 

We can calculate the single photoelectron rate for CT1 and MAGIC. Note 
that the poor angular resolution (a single pixel with FOV of 0.25° for CT1 
and 0.1° for MAGIC), is compensated for by the large mirror area A (10 m 2 
for CT1, 250 m 2 for MAGIC). 
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The expected photoelectron rates are calculated with the following expression: 
R = J Fx-vQu- A ■ (3 du, (1) 



where 77 is the mirror reflectivity and Q v the frequency dependent quantum 
efficiency of the PMs. The parameter j3 takes into account the photoelectron 
detection on the first dynode (0.9) and the light guide efficiency (90%). 



For the Crab pulsar, the theoretically expected rates (Eq. 1) are given in Table 
2, for CTl and MAGIC telescopes. 



Telescope 


Total Pulsar Rate 
(phe -s~ 1 ) 


Nebula Rate 
(phe -s" 1 ) 


LONS Rate 
(phe -s -1 ) 


Background Rate 
(phe -s" 1 ) 


CTl 


1.5 x 10 5 


1.4 X 10 8 


3.3 x 10 7 


1.7 X 10 8 


MAGIC 


3.7 x 10 6 


3.5 X 10 9 


1.3 x 10 s 


3.6 X 10 9 


Table 2 



Expected currents (in photoelectrons/s) of the optical emission of the Crab pulsar 
in the central pixel of the CTl and MAGIC telescopes. Background rate means the 
total background: LONS + Nebula 



If we assume that the PM has no inherent fluctuations (for example the fluctu- 
ations of the electron multiplication process between dynodes), we may derive 
the basic scaling parameter X = p^J{N) based on photon (Poissonian) statis- 
tics, for all test statistics, which measures the integrated difference between 

n 

the measured pulse profile and a uniform distribution (no pulsed signal) (J3(). 
The parameter p is the number of pulsed events, divided by the total number 
of events N (signal plus noise). This X also corresponds approximately to the 
expected DC excess (measured in Gaussian standard deviations). Assuming 
that the ephemeris for folding the arrival times of optical photons from the 
Crab pulsar to phases is correct, the expected scaling parameters for CTl and 
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MAGIC (taking Ilons from (|7|), although we expect it to be higher), based 
on photon statistics alone, are: 



From the basic scaling parameter we infer the expected time to detect the 
optical pulsation for a given significance. The theoretical rate from the Crab 
pulsar shows that it would be possible to detect the Crab optical pulse with 
the CT1 telescope in T=75 ms, with a significance of 3a. For MAGIC, the 
situation is even better because the main pulse could be detected within 5 ms, 
which is the approximate duration of this pulse. This implies that MAGIC 
has the capability to monitor single pulses from Crab, based on theoretical 
arguments, and using the LONS intensity corresponding to an extragalactic 
source. 

This result is only valid if the photodetector itself has no noise, and the only 
source of fluctuations is Poissonian in nature, due to the signal. The expected 
detection times increase dramatically when we take fluctuations into account, 
as shown by the measurements presented below. 

3 Measurements 

Measurements were done with the HEGRA CT1 telescope, during February 
ll-16th 2002 and during November l-8th 2002. The standalone HEGRA CT1 
telescope has a reflector area of 10 m 2 and a camera consisting of 127 PMs, 




(2) 




(3) 



each with 0.25° FOV. 
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Datasets were obtained with the Crab in the central pixel, and OFF source 
runs were also taken to check for any systematics. For optical observations the 
central pixel PM was modified. The DC branch, designed to monitor the DC 
current of the pixel, was adjusted to detect pulses of ~3 ms (timescale of the 
pulse width), and the AC branch, which is designed to transmit the ns fast 

Jl 

signals generated by the Cerenkov showers was removed (J5|). 

The DC current was measured using a National Instruments PCI NI 6043E 
digitising card. The sampling rate was set to 2 kHz and the pixel line was AC- 
coupled to reject low frequencies, and increase the card's effective dynamic 
range. 

4 Analysis and Results of the Pulsed Optical Signal Measurements 

The arrival times £j were folded to obtain the phases $j. 

$i = $ + v(u - 1 ) + \y{h - to) 2 , (4) 

where is, is are the frequency and frequency derivative, and to is the reference 
epoch. 

To transform these arrival times to the solar barycentre system, we used the 
TEMPO software Jf]). A correction due to a small average drift in the digi- 
tisation card clock, as measured with respect to an atomic reference clock, 
was also introduced. The Jodrell Bank ephemeris (frequency, first frequency 
derivative and reference time) applicable to these observations are given in 
Table 3 (Jodrell Bank Crab pulsar monthly ephemeris). 
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to 



V 



V 



Date 



(Hz) 



(MJD) 



(1Q- 15 



sec 



15th October 2002 



52562 



29.8132421622 



-373777.06 



15th November 2002 



52593 



29.8122410972 



-373758.90 



Tab: 



e 3 



Radio Timing Parameters for the Crab pulsar for the indicated epochs. 

A phaseogram was produced for each independent frequency confined to the 
wide frequency range between 26.0 and 32.0 Hz. The independent Fourier 
spacing (IFS) between this range of independent frequencies is known to be 
1/T (periodogram bin size), where T hours is the observation time (2 hours). 

The folded intensities for each test frequency were tested against a uniform 
distribution by performing a x 2 fit to a constant intensity. The reduced x 2 
was calculated (Figure 1). A maximum value of the reduced x 2 °f ~ 154 
was found and the periodogram signal is contained within one IFS of the 
expected Jodrell Bank ephemeris, consistent with the uncertainty in the drift 
expected from the clock on the digitisation card. Absolute phase matching with 
the expected radio/optical/X-ray/7-ray main pulse position was not possible, 
since an absolute UTC match was not done, and the drift was not accurately 
measured during observations. 

Furthermore, the frequency resolution for the 2 hour observation was too 
coarse to resolve the differences between the ephemerides corresponding to 
15 October and 15 November. The predicted frequency shown in Figure 1 is 
therefore consistent with both ephemerides. 

Figure 1 shows our measured Crab pulse profile, with two sharp peaks, with 
a ~ 0.42 phase difference. The corresponding pulsed fraction can now be used 
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to isolate the LONS and cross-check the previous measurement (|7|). Using our 
results we can infer the MAGIC response for the Crab pulsar. 
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Fig. 1. Observations of the optical Crab pulsar with CT1. On the top, the two peaks 
are clearly visible in the light curve with a separation in phase of ~ 0.4- On the 
bottom (b) the excess in the reduced x 2 searching in a wide range of frequencies, using 
the 15th October ephemeris. The difference between the frequency found with 15th 
October ephemeris and 15th November is consistent with the drift in the digitization 
card clock. The inset shows the value of the reduced \ 2 vs. x = (v — v ) / IFS , where 
v Q corresponds to 15th October. 
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5 Determination of the galactic anticenter LONS 

In this section, we will determine the galactic night sky background around 
the Crab pulsar by using the pulsed fraction of the signal detected in CT1. 
This result will be checked by comparing it with the LONS calculated from 
the intensity of the CT1 camera response to ( Tauri. 

It is possible to derive the pulsed fraction of the optical Crab signal from 
the density function of the light curve. This density function results from the 
unpulsed background (nebula and night sky background) with fraction 1 — p 
and the density function f s {9) of the pulsed signal coming from the Crab 
pulsar: 

f(0)=pfs(e) + (i- P ), (5) 

where p is the pulsed fraction, if the density function is normalised to unity. 
The off-pulse interval (region of minimum intensity) following the interpulse 
and preceding the main pulse was used to estimate the pulsed fraction. 

Taking the proper errors in each bin into account (exploiting the Central Limit 
Theorem), we derive a pulsed fraction for the Crab of p — (5.00 ± 0.02) x 10~ 4 . 

To check the stability of the signal with time, the total data set of two hours 
was divided into six data sets of 20 minutes each. The pulsed fraction, as well 
as the ratio between the first and second peaks was calculated for each subset. 
The results are shown in the Figure 2. It is clear that both parameters appear 
to be stable within statistical errors. An indirect indication of the stability 
of telescope tracking during two hours observations is shown by the constant 
value of the pulsed fraction. 
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Pulsed fraction * 10 4 



— t — — + — 
Ratio (2 peak/ 1 peak) 







5 6 7 
20 min sets 



Fig. 2. Test of the telescope response to the pulsed signal in two hours. The pulsed 
fraction times 10 for the set of 20 min observation slices is at the top of the plot. 
The ratio between the first and second peak remains constant with time, as shown 
at the bottom. 

There is a difference between the measured CT1 Crab pulsed fraction of 5 x 



10 4 and the theoretical value given in the previous section (p 



theor 



8.6 x 10" 4 1 



The discrepancy can be attributed to the uncertainty on the measurement of 
the LONS, the influence of the mirror point spread function (PSF), and the 
fact that the Crab pulsar was off-center due to uncorrected but constant mis- 
spointing in the telescope. 



To determine the galactic LONS around the Crab pulsar, we first need to 
know the PSF of the telescope. With this purpose, we used the bright source 
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( Tauri to measure the amount of light, which falls within the central pixel. 
The same source is used to measure the telescope misspointing. 

The selected source was ( Tauri (ICRS 2000.0 coordinates are 5 h :?>7':?>&.7§" 
RA, 21°:08':33.0" DEC with a magnitude of 2.84), which also appears in the 
camera field of view when CT1 points to the Crab pulsar. The advantages 
of using this reference star are the following: (1) we know its photometric 
properties and (2) the atmospheric extinction will be the same for both the 
pulsar and the star, given their proximity in the sky. The PSF of the telescope 
was determined from a raster scan with a resolution finer than the 0.25° pixel 
size (see Fig. 3(a)), with zero position corresponding to the true position of the 
star). Since the telescope is focused to ~6 km (maximum Cerenkov shower), 
the image of a point-like source will be smeared out around 0.14°. 

The signal from ( Tauri is slightly shifted to lower right ascension values, 
due to the telescope mispointing. This was accounted for by calculating the 
maximum value of the signal, and measuring the difference between the shifted 
measurement and the true position. A schematic figure of the CT1 central pixel 
is also drawn in Figure 3. To calculate the true signal collected by the central 
pixel, the signal intensity was fitted to a bidimensional Gaussian (fig. 3 (b, 
c)), and the fraction of the signal within the central pixel was measured to be: 



From the misspointing correction factor II, it was possible to derive the LONS 
rate from: 



IT = 



■r pixe i j2k q ^ ^ drde 

JS°jS"G(r,9)drdB 



= 69%. 



(6) 



P = 



■pulsed 



(7) 



Rpulsed + Rnebula + RlONS ' 
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b) 




Fig. 3. DC currents due to Q Tauri on the CT1 central pixel, in 0.1 fiA units. Left 

(a) : C Tauri signal in the CT1 camera, with background already subtracted. The 
zero position corresponds to the true position of the star in equatorial coordinates 
(maximum intensity for zero mispointing). The approximate shape of the central 
pixel is drawn to illustrate the size of the PSF relative to the pixel size. Top right: 

(b) C Tauri observed intensity in a 2-dimensional plot. Bottom right: (c) The same 
signal, fitted with a bidimensional Gaussian. 

where p is the empirical pulsed fraction, R pu ised is the expected rate (table 2) 
of Crab pulsed events, R ne buia the expected rate from the Crab nebula (table 2) 
and Rlons is the galactic background rate around the Crab nebula. 

Finally, assuming that the contribution of the pulsed signal is negligible com- 
pared to the background due to the nebula and the LONS we have: 

= J F^ aed ■ V ■ Q ■ A^e^Udu 

P j pnebuia . r] . q . A ctl e yv Udu + J F^ns . Q . rj . Q . A dl e lv du ' 1 j 

where rj is the mirror reflectivity, Q the quantum efhcency of the PMs, A ct \ is 
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the CT1 mirror area, II is the fraction of light from a point source captured by 
the central pixel, Q is the solid angle subtended by the central pixel, and e 7 „ is 
a conversion factor between number of photoelectrons and voltage measured 
by the digitization card, which cancels out in the above expression. 

Therefore, using the pulsed fraction of the total Crab signal, the frequency in- 
tegrated LONS (weighted with the frequency dependency of the PM quantum 
efficiency) is: 

Ilons = J F^ ONS dv = (4.6 ± 0.1) x 10 12 ph ■ m- 2 ■ s' 1 ■ sr' 1 (9) 

To cross-check this result, we used the signal detected by CT1 from ( Tauri in 
the central pixel (fig. 3): If i cp is the total observed current intensity from the 
direction of ( Tauri in the central pixel, Tauri is the expected contribution 
to the central pixel current from ( Tauri, and e 7 « is the conversion factor from 
photoelectrons to current (in //A), which is calculated from 

e 7i n J F^ Taun -rj-Q ■ A ctl dv = i c Tauri = i cp - i LON s ~ i cp - (10) 

The last approximation is valid since ( Tauri is a very bright star rendering 
the contribution from the LONS in the central pixel negligible. F^ Tauri is 
tabulated including atmospheric extinction. We obtain a measured value of 
e 7 j = 2.95 x 10~ 9 mA ■ s ■ phe' 1 . From this value we obtain the LONS by 
subtracting the PSF of ( Tauri from a square degree image centered on this 
star (shown in Figure 3). The error on the LONS is minimal near the outer 
wings of the PSF, since the LONS is dominant in that area. An annular region 
is used to estimate the LONS, giving 

Ilons = J F^ ONS dv rsj (4.3 ± 0.1) x 10 12 ph ■ m~ 2 ■ s" 1 • sr' 1 , (11) 
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after normalising to the solid angle subtended by the selected annular region. 



The two values for the galactic Ilons ( e Q- 9 & 11) calculated by two in- 
dependent methods agree within ~ 7% and are ~ 2.7 times higher than the 
previous measurement J?|), used for our theoretical predictions as we expected. 
The LONS intensity calculated from the pulsed fraction is the most accurate, 
being the correction factor II uncertainty the highest, since the LONS deter- 
mination from £ Tauri also includes the uncertainty in the background region 
around the star. 

6 Sensitivity for Pulsed Detection in Optical 

The sensitivity of CT1 telescope to detect the Crab optical pulse can be ex- 
pressed by the " standard deviation" s of the signal relative to the behavior of 
the reduced x 2 a t frequencies different from the Crab frequency: 



for each observation time interval, where ^background is the mean value of x 2 red 
at such offset frequencies and (Tbackground the corresponding standard deviation, 
which is constant for all selected times T. The total 2 hours observation time 
was split into L = 7200 s/T independent datasets, each of length T, from 
which the average value of the statistic s was calculated for T = 50 s to 7200 s 
in steps of 50 s. It can be shown that this average value of s is: 



s = 




(12) 



^background 



s ~p 2 ■ T ■ g(k, f s ) 



(13) 
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where k is the number of bins in the phaseogram, f s the signal pulse shape 
and g(k,/ s ) a function depending on these quantities. Note that s is not a 
Gaussian standard deviation, since the number k is too small, but we expect 
that < s >oc a 2 , where a ~ p • VT would be the Gaussian significance of a 
DC detection of the same signal. Figure 4 shows a graph for log < s > vs. 
logT, showing a linear relationship between s and T as predicted. 

Using the calculated value of the galactic LONS from the pulse analysis, we 
estimate the required time for MAGIC to detect the Crab pulsar in the optical, 
given the much larger mirror area, improved quantum efficiency, and pixel size 
(0.1°). By scaling the graph of log< s > vs. logT accordingly, Figure 4 shows 
that MAGIC can detect the Crab pulsar in 30 s at the 5a level. The line on 
top shows the prediction for a PM without inherent fluctuations (eq. 3), using 
the galactic LONS calculated. 

7 Conclusions 

The Crab pulsar was measured in the optical with the HEGRA CT1 telescope 
(the smallest Cerenkov telescope making this detection). A 5<x detection re- 
quires an observation time of 200 seconds. 

Since the design of the CT1 PM base is very similar to that of MAGIC, we 
extrapolate CT1 results and get a 5<r detection with MAGIC (~ 25 times larger 
dish area than CT1) in a substantially shorter time of 30 sec. Considering that 
it is possible to reduce the LONS using an array of HPD's, combined with the 
superior PSF of MAGIC, it could still be possible to resolve single pulses. 

The use of a PM base with the technical modifications described above is 
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log(time(s)) 



Fig. 4. log of "standard deviation" s of the signal related to the background versus 
observation time. The CT1 empirical data are represented with filled bullets while 
the empty ones represent the extrapolation for the MAGIC telescope. The horizontal 
line shows the "5a" level. 

incompatible with normal VHE observations. The modification could only be 
considered for the central pixel of the telescope camera. 

The galactic LONS background on the surroundings of the Crab region was 
calculated using the Crab pulsed fraction. Only 7% disagreement was found 
with a second method using a bright reference star and a different background 
region (still within the camera field-of-view) . Our results are however on aver- 

Jl 

age 2.7 times higher than the extragalactic LONS from ((3). The advantage of 
this new method is that it relies on the stability of the clock, and is indepen- 
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dent of other sky-techniques in the DC mode, provided that the PSF is known. 
The determination of the LONS for the Crab region is important, since the 
threshold energy determination and the spectrum depend on the level of the 
LONS. 

Simultaneous Crab pulsar optical data and Cerenkov VHE observations could 
be very important for monitoring atmospheric extinction. In fact, high qual- 
ity observations should show a very specific zenith angle dependency of the 
pulsed fraction. Thus, the "seeing condition" at a given time may also be 
derived from the observed pulsed fraction corrected for atmospheric extinc- 
tion. Furthermore, apart from using optical observations to provide a timing 
ephemeris for VHE observations. 

We can also monitor deviations from the normal cubic law for pulsar spindown 
(assuming the pulsar braking index is detectable) as a result of timing noise 
and glitches. Thus, frequent monitoring of the Crab by MAGIC (for calibration 
purposes) should also improve the chance of detecting glitches. In fact, very 
specific software can be written to test glitching behavior during observations. 
A specially designed central pixel detector with high quantum efficiency, and 
relatively low noise will be very useful for monitoring the state of optical 
emission from 7-ray sources such as cataclysmic variables, X-ray binaries, 
microquasars, AGNs and possibly GRBs. 
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